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ABSTRACT 


A  FORTRAN  program  is  explicitly  listed  for  the  basic  integral 
solution  of  the  wake  of  a  hypersonic  blunt  or  conical  reentry  body. 
The  near  wake  pressure  gradient,  nonequilibrium  chemistry,  and 
laminar -turbulent  transition  are  included. 
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SYMBOLS 


Text  Program* 

a 

A', A" 


CD 


ALT 

CD 


CP 


D 

^ijave 

d  LD 

E 


^chem 


eint 


g 

H 


*See  footnote  at  end  of  symbols. 


Definition 

Base  area,  feet2,  Equation  (70). 

Constants,  Equations  (60) 
and  (61). 

Altitude,  feet 

Drag  coefficient. 

Equation  (70). 

Specific  heat  at  constant 
pressure  of  air  mixture, 

“K,  Equation  (15). 

Specific  heat  of  species  i, 

°K,  Equation  (15). 

Drag,  pounds,  Equation  (69). 

Average  binary  diffusion 
coefficient  between  species 
i  and  j.  Equation  (8). 

Base  diameter  of  body,  feet. 

Electronic  energy,  feet2  per 
second2.  Equation  (54). 

Zero-point  energy,  feet2  per 
second2.  Equation  (54). 

Internal  energy,  feet2  per 
second2  per  unit  mass. 
Equation  (4). 

Vibrational  energy,  feet2  per 
second2  (equal  to  Planck's 
constant  times  y  <«). 

Electronic  degeneracy. 
Equation  (54). 

Total  enthalpy,  feet2  per 
second2.  Equation  (57). 


IV 


Text 

Program* 

Definition 

h 

LHXN  (h(x,n)) 

Enthalpy,  feet2  per  second2. 
Equation  (53). 

Hnn0 

HNNO 

O2H/9n2)0  ,  Equation  (43). 

»v 

I 

Unit  tensor,  Equation  (6). 

K 

Thermal  conductivity,  - 
Equation  (7). 

k 

Boltzmann's  constant. 

Equation  (54). 

Le 

LEO  (Le0 ) 

Lewis  number,  Equation  (10). 

M 

W 

Molecular  weight  of  atmo¬ 
sphere  model.  Equation  (74). 

Moo 

MIN 

Free-stream  Mach  number, 
Equation  (64). 

m 

Transformation  variable, 
Equ-^ion  (23). 

M(i) 

FM 

Mass  of  i*h  species. 

Equation  (54). 

n 

Transformed  radial 
coordinate.  Equation  (24). 

P 

PRESS 

Pressure,  pounds  per  foot2, 
Equation  (6). 

Pr 

PRO  (Pr0 ) 

P r and tl  number.  Equation  (10). 

q 

Dynamic  pressure,  pounds  per 
foot2.  Equation  (70). 

q 

Heat  flux  vector.  Equation  (4). 

R 

Radial  distance  to  shock,  feet. 
Equation  (69). 

r 

Radial  coordinate,  feet. 
Equation  (25). 

T 

TEMP 

Temperature  (kinetic),  °K 
Equation  (75). 

Tm 

TEMPM 

Molecular -scale  temperature, 
Equation  (76). 

❖See  footnote  at  end  of  symbols. 

v 


Text 


Definition 


Procram* 


START  X 
STOP  X 


ZOE, 


Time,  seconds. 

Velocity  in  downs tream  (x) 
direction,  feet  per  second. 
Equation  (9). 

Free-stream  velocity,  feet 
per  second. 

Diffusion  velocity,  feet  per 
second.  Equation  (5). 

Velocity  in  the  radial  (r) 
direction,  feet  per  second. 
Equation  {9). 

General  fluid  velocity,  feet 
per  second.vector, 

Equation  (1). 

Downstream  coordinate,  feet. 
Equation  (9). 

Laminar  turbulent  transition 
coordinate,  feet.  Equation  (64). 

Increment  of  x,  feet. 

Initial  x  value,  feet. 

Terminal  x  value,  feet. 

End-point  altitudes,  feet. 
Equation  (76). 

Transition  coordinate  para¬ 
meter,  Equation  (64). 

End-point  temperatures. 
Equation  (77). 

Compressibility  factor  on  the 


1  ;M 

mm 


'  j'S. 


- 


A,  6,  C,  D,  E, 
F.K,  L,  N,  P 

B,  C,  D,  Ej,  Hj 


ALPH 


*See  footnote  at  end  of  symbols. 


Constants  (Paragraph  4), 

Atmosphere  constants. 
Equations  (77)  through  (81). 

Species  concentration. 
Equation  (11). 


Text 

Program* 

Definition 

6 

DEL 

Radial  distance  in  x,  r 
coordinate  system,  feet, 
Equation  (72). 

6m 

DELM2  (6m) 

Transformed  distance  to  outer 

edge  of  wake,  feet,  Equation 
(25). 

0 

THETA 

Momentum  thickness,  feet. 
Equation  (17). 

eE 

Energy  thickness,  feet, 

Equation  (18). 

8(i> 

Concentration  thickness, 
feet,  Equation  (19). 

^C 

DELTAC 

Constant  depending  on  initial 
conditions,  Equation  (44). 

A  W 
e 

Defined  by  Equation  (59). 

P 

MU 

Viscosity  coefficient,  — — z~-—  , 

V(i) 

Equation  (6). 

Vibrational  frequency  of  i**1 
molecular  species,  Equation 
(54). 

p 

DENSG 

lb-sec2 

Density,  Equation  (1), 

4>n 

PHI 

Cone  half-angle,  degrees. 

4> 

S2B 

Shock  half -angle,  degrees. 
Equation  (67). 

4* 

*  *  * 

Coordinate  in  initial  cylindrical 
coordinate  system  (not  used 
in  equations). 

i 


c 


SUBSCRIPTS  AND  SUPERSCRIPTS 

I  Index  for  eight  chemical 

species. 

C  Initial  x  coordinate  for 

calculation. 


*See  footnote  at  end  of  symbols. 
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vii 
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Text 

e 

0 

oo 


Program* 

E 

0 

IN 


Definition 

Values  at  n  =  1 . 

Values  at  n  =  0, 

Free -stream  conditions. 


! 


u 


r 

r 


if 

r  f  ^ 

Jf 


*Only  a  few  representative  program  symbols  are  listed;  however, 
unlisted  symbols  can  be  deduced  from  the  text  equations.  The  units  used 
in  the  program  are  feet,  pounds,  seconds,  and  degrees  Kelvin.  Corres¬ 
ponding  units  are  assumed  throughout  the  text  unless  noted  otherwise. 


vm 
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1 .  Introduction 


The  analysis  presented  in  this  report  is  essentially  a  program¬ 
ming  of  the  hypersonic  axisymmetric  wake  problem  as  formulated  by 
Bloom  and  Steiger  1  and  Ness  and  Fanucci. 2  The  object  of  this  v/ork  is 
to  use  the  fundamental  wake  analysis  as  a  basis  for  an  extension  of  the 
problem  to  include  contributions  from  ablation  and  base  mixing  and  to 
determine  their  effects  on  the  initial  wake  conditions  and  concentrations 
of  radiating  species.  The  FORTRAN  program  is  in  Appendix  A.  To 
be  complete  and  internally  consistent,  the  program  is  described  in 
detail  beginning  with  the  basic  conservation  e^-’ations.  The  solutions 
are  obtained  for  an  air  model  of  eight  species  ■„  ->ing  nonequilibrium 
chemistry.  The  downstream  pressure  gradient  in  the  wake  and  the 
laminar -turbulent  transition  are  included.  However,  a  number  of 
important  assumptions  have  been  employed  to  make  the  analysis  tenable. 
The  units  used  in  the  program  are  feet,  pounds,  seconds,  and  degrees 
Kelvin.  Corresponding  units  are  assumed  throughout  the  text,  unless 
noted  otherwise. 


2.  Definition  of  the  Flow  Field 


The  basic  regions  of  the  flow  field  accompanying  a  hypersonic 
body  as  it  enters  the  atmosphere  have  been  extensively  studied  both 
experimentally  and  theoretically.  Air  is  heated  irreversibly  and 
compressed  by  the  strong  detached  bow  shock  associated  with  a  blunt 
body.  The  fluid  entering  the  boundary  layer  overexpands  around  the 
rear  edge  of  the  body.  In  this  base-flow  region,  the  fluid  is  recom¬ 
pressed  and  turned  back  parallel  to  the  axial  direction  producing  a 
secondary  shock  which  is  weaker  than  the  main  bow  shock.  Fluid  which 
is  not  captured  by  the  recirculating  base  flow  passes  through  a  high 
pressure,  minimum  diameter  neck  and  into  the  viscous  wake.  A  conical 
body  is  characterized  by  an  attached  shock  which  produces  much  less 
heating  in  the  inviscid  portion  of  the  wake  and  by  a  thicker  boundary 
layer. 

The  inner  viscous  wake  may  be  divided  into  two  regions  distinguished 
by  the  presence  of  an  axial  pressure  gradient.  The  near  wake  extends 
from  the  vicinity  of  the  neck  to  the  downstream  coordinate  at  which  the 
pressure  has  essentially  decayed  to  the  ambient  value.  The  far  wake 
is  then  controlled  primarily  by  diffusion  processes.  Large  pressure 
gradients  and  high  temperatures  characterize  the  initial  portions  of  the 
near  wake.  The  outer  inviscid  wake  passes  through  the  bow  shock  and 
the  recompression  shock  to  be  "swallowed"  by  the  radial  growth  of  the 
inner  wake  at  some  downstream  coordinate  determined  by  the  flight 


i 


characteristics  and  the  laminar -turbulent  transition  point  of  the 
inner  wake. 


3.  Basic  Conservation  Equations,  Transformations,  and  Boundary 
Conditions 

When  both  the  chemical  and  fluid -mechanical  properties  of 
the  flow  are  considered,  conservation  relations  can  be  written  for  the 
momentum,  energy,  and  chemical  species.  In  addition,  an  equation 
expressing  global  continuity  (conservation  of  mass)  is  valid.  Using 
concise  vector  notation,  the  continuity  equation  is 

ff  +  V  •  (pV)  =  0.  (1) 

Because  the  calculations  are  done  under  the  assumption  of  steady  flow, 
all  time  derivatives  vanish. 

V  •  (pv)  =  0  (2) 

Similarly,  the  momentum,  energy,  and  species  equations  are 


respectively 

V  •  (pW  +  p)  =  0 

(3) 

V  •  ({p|vpv  +  peint  V  +  V  •  p  +  q)  =  0 

(4) 

V  •  p(i)(v  +  V(i))  =  fel-)  , 

p  \  3t  /chem 

(5) 

p  is  the  pressure  tensor,  is  the  heat  flux  vector,  and  is  the 
diffusion  velocity. 

p  =  pi  -  p  V v  +  (Vv)t  -  (v  *  v)l]  (6) 

q  =  -KVT  +£  pWvWhW  (7) 

i 


2 


V<i>  =  -Dij 


Va(i)/a(i) 


(8) 


ave 


Superscript  i  represents  the  chemical  species  of  the  reacting  air 
mixture.  ^ijave  the  averag®  binary  diffusion  coefficient  between 

species  i  and  j.  represents  the  unit  tensor,  K  is  the  thermal  con¬ 
ductivity,  and  eint  represents  the  internal  energy  per  unit  mass  of  the 
air  mixture.  The  remaining  notation  can  be  obtained  from  the  list  of 
symbols. 

The  physical  coordinate  system  consists  of  a  downstream  coordi¬ 
nate  x,  a  radial  coordinate  r,  and  an  angular  coordinate  ({>.  However, 
the  cylindrical  symmetry  of  the  problem  removes  any  <}>  dependence. 

It  is  also  possible  to  apply  boundary  layer  approximations  if  the 
gradients  normal  to  the  flow  direction  are  stronger  than  those  along 
the  flow  direction. 


a.  Momentum 


h. +  £(ri>uv>  -  £ 


I  dAj.  9 

\rfiar/+ax 


(rp)  =  0 


(9) 


b.  Energy 


_a_ 

3x 


(rpuH)+A(rpvH)=^(^r|f) 


9r 


£<Pr-»ru£ 


+  £[£(Le-1)r£ h<i>  If 


(10) 


(i) 


c.  Chemical  Species 


^(rpuaW)  +  ~(rpva(l)) 

3  /  Le  3a(i)\  /  3p{i)\ 

~3r  v'Pr  3r  J*  VdT  )chem 


(ID 


1 


the  integral  form  of  the  momentum,  energy,  and  chemical  species 
equations  in  the  x,  r  system  is 


^(lnPeu|)  -  (J  ^  rdr^  in  u( 


-Si-  |E.  0 

2peue  dx 


Peue0E  =  constant 


|i(l>  +  e'1'  A. 

dx  dx 


<£n  *  Uel  -  (  /  rdr) 


—  s{^—] 

;Ue  *  \  at  /chem 


The  transformed  radial  coordinate  n  is  now  introduced  by  the  Dorodnitzen 
transformation  through  the  relationships 

mdm  =-£-  rdr  (23) 

Pe 

m  =  n6m  (24) 

which  lead  to 


6  2  ndn  =  —  rdr.  (25) 

m  pe 

The  effect  of  this  transformation  is  to  remove  the  radial  density  depen¬ 
dence  from  the  calculations  by  introducing  the  x,  n  coordinate  system 
in  place  of  the  cylindrical  x,  r  system. 

Several  sets  of  boundary  conditions  apply  to  the  concentrations, 
velocity,  and  enthalpy  at  the  axis  (n  =  0  )  and  on  the  outer  edge  (n  =  1). 


I  •:& 


H  =  H‘:  (ir)e*  (0)e'° 


These  boundary  conditions  follow  from  a  consideration  of  the 
physical  model  of  the  symmetrical  wake.  From  these  conditions  it  is 
possible  to  formulate  profile  expressions  of  the  proper  form 


a  ^  (x,  n)  -  a0  ^  (x) 

(i)  (i)  ,  , 

a  (x)  -  a0  (x) 
e 


u(x,  n)  -  u0  (x) 
ue  (x)  -  u0(x) 


6n2  -  8n3  +  3n4. 


(32) 


The  momentum  Equation  (9),  the  energy  Equation  (10),  and  the  chemical 
species  Equation  (11)  are  next  specialized  to  n  =  0  and  n  =  1  after 
transformation  to  the  x,n  system. 


6 


7 


(41) 


Pe  ue^E  =  (Pe  ue®E)c 
where 


The  shape  parameter  Hnn()  is  then 


(42) 


Hnn0 


4(He 


H0 ) 


349  ue  +  311  u0 
43  ue  +  23  u0 


(Pe  Sm)  c  Ac 
Pe  6m  (43  ue  +  23  u0) 


(43) 


Ac  depends  on  the  initial  conditions  as  indicated  by  the  subscript  c. 

Ac  =  4  (He  -  H0c)  (349  ue  +  311  u0)c  -  (Hn„0)c  (43  ue  +  23  u0  )c  (44) 


4.  Reacting  Air  Mixture 

A  reacting  air  mixture  containing  eight  chemical  species  is 
used  to  represent  the  fluid  chemistry.  The  air  species  are  N2 ,  02, 
N,  O,  NO+,  02  ",  NO,  and  e” .  The  22  allowed  chemical  reactions 
involving  these  eight  species  are 


(R1  -  R5) 

Oz  +  X  ~  20+  X 

(X  =  O,  02  ,  N,  N2  ,  NO) 

(R6  -  RIO) 

Nz  +  X  Z1  2N  +  X 

(X  =  O,  02  ,  N,  N2  ,  NO) 

(Rll  -  R15) 

NO+xr  N+O  +  X 

(X  =  O,  02  ,  N,  N2  ,  NO) 

(R16) 

NO  +  O  ^  N+  Oz 

(R17) 

N2+  O  ~  NO+  N 

(R18) 

N  +  O  ^  NO+  +  e- 

(R19) 

NO+  e-.^T  NO  +  +  2e- 

(R20) 

Nz+  02~  2NO 

(R21) 

02  +  e  "  +  02  ■*—  Of  +  02 

(R22) 

NO  +  +  Oz  -  ^  NO  +  02 

Constants  involved  in  the  chemical  kinetics  are  given  by  Ness  and 
Fanucci.  2 


The  chemical  kinetics  enter  the  analvsis  through  the  chemical 


production  terms 


(Hr-) 


chem 


in  the  law  of  mass  action.  The  partition 


8 


1 


functions  are  derived  by  assuming  equilibrium  for  all  degrees  of  freedom. 
The  resulting  expressions  are 


1  1*0™]  .  .  m<°> _ m<°> 

pHWcW  '  WC-D-E>-JNq 


2 

P  \  3t  /ch 
i  ( a  p 

p  \  at  /c 


(N)  (N) 

__  m  .  __  .  m 

=  B  +  "  (NO)  (C  *  D  +  E)  (NO+)  F 
chem  m  rn 


(02)  (02  ) 

m  ,  _  ,  _  ,  m 


-  (N-pi>  <47> 

P  \  at  /chem  (NO)  2  (°2 ) 


i  /a  (N*M  m(N2) 

7  \~Vr~ /chem  =  "B  “  (E  +  *L) 


1  9p 


(NO)' 


(NO) 


—  \  ZJZ - J  =E  +  L-  C-  D-  K  +  P1 

p  \  3t  /chem  (Oz  ^ 


1  (  a 


?). 


(NO+ )  (NO+ ) 


- I  =  F  +  - -  K  -  — _  P' 

at  /chem  (NO)  (02) 


1  (  9P  )  _  m 

p  \  at  ) chem  (NO+ ) 


(NO)  _(02) 


l( L^.) 

’P  \  9t  /chem 


=  N  -  P‘  . 


The  constants  A,  B,  C,  D,  E,  F,  K,  L,  N,  and  P'  are  given  by  Ness 
and  Fanucci.2 

The  total  static  enthalpy  h  of  the  air  mixture  is 


-Za(i>  h(i>  .  w-S? 


The  individual  species  states  enthalpy  is  a  function  of  the  static  tempea 
ature  and  molecular  and  atomic  constants  only. 


-  -if- -f .. 


9 


i 

£ 

%- 


£ 


K 

*r. 


f-T 

r- 


L‘ 

£ 


Ji)  1  ,kT_  g  (l) 

"2  (i)  (i)/  _  (i) 

m  m  I  Evit> 


-fvib.  .  - 


m 


(i) 


-Ejg(i) 

S  g^(l)  E^(i)  e  ^ 


-E 


(i) 


Sgi 

n 

From  the  equation  of  state, 

T“>  M 

p  =  p  kTE-V 

•  m 


(i)  P  kT 


E  chem 

(i) 


m 


(54) 


(55) 


(56) 


i 


.$piggj 


n  wq 

\-f  z~' 

'  *"* 


'-¥£3~jp^za 


The  temperature  derivative  can  be  obtained  from 

He  =£aeU)  h=(i>  +  4  <57) 

i 

Substitution  of  the  he^  values  into  Equation  (57)  and  differentiation  gives 

A” 


d£n  T 


A'  + 


e_ 


kT, 


dx 


V*  ae(i)  \  U  {i))  2  I  a  (i)  Evib(i)  7]  °e(l) 

2Jj?ri[e  1  e  m<« 


(i) 


Eyjb(l) 

kTe 


Evib(l) 

kTe 


(58) 

(59) 


-1 


10 


(60) 


dae(i),  5_  1  dge(i) 

dx  ZmJ  2  (IT  dx 


A" 


Severn 


i 


(61) 


5.  Laminar -Turbulent  Transition  and  the  Pressure  Gradient 

The  laminar -turbulent  transition  is  treated  by  changing  the 
viscosity  model.  The  laminar  portion  of  the  wake  is  governed  by  the 
Sutherland  equation  . 


-8  T  $IZ 

p0  =  3.05  x  10 

T0  +  110 

The  turbulent  viscosity  is 

po  =  0.02  6m  pe  (ue  -  u0 ) 


(62) 


(63) 


The  transition  coordinate  is  computed  using  the  equation  formulated 
by  Zeiberg.3 


xt-ii*JL 

Poo  U°° 

Y  is  defined  by 


5  6 

+  — 


sin2  <j> 


M2,  sin2  <J>  +  5  \7M*  sin2  <J>  -1 


Moc  <  8  :  log10  Y  =  4.  744  +  0.  124MQC  +  0.  00976  M<g 
Moc>  8  :  log jq  Y  =  5.  48  +  0.  1 1  Moc 
A  cubic  can  be  written  for  the  shock  angle  $ 

(sin2<J>)3  ^  +  1.4  sin2  4>cj  (sin2  4>fj)2 

i+  (1-44+l^)sinJ^sinI*-;i^’s0 


ZMl 


MJo 

is  the  cone  half -angle 


r 


(64) 

(65) 

(66) 


(67) 


11 


/ 


r*s\ 


vat 


6.  Sequence  of  Calculations 


A  printout  of  the  FORTRAN  source  program  is  given  in  Appendix 
A.  In  the  right-hand  margin  references  are  made  to  the  equations  used 
in  the  calculations.  Important  program  symbols  are  defined  in  the  list 
of  symbols  and  correlated  with  the  symbols  used  in  the  text. 


The  input  data  include  a  values  for  02  ,  N2  ,  O,  N,  NO,  02~,  NO"*', 
and  e“  and  ae  values  for  02  and  N2  .  The  remaining  initial  ae  values 
have  been  set  equal  to  zero  internally.  Initial  u0 ,  H0  ,  T0  ,  4>n,  Z0  , 
Hnn0  >  and  x  values  as  well  as  Pr0  ,  Le0  ,  Cq,  d,  ,  and  final  x  values 
are  used  as  input.  The  altitude  must  be  provided  to  initiate  the  atmo¬ 
spheric  subroutine  for  the  calculation  of  free -stream  temperature, 
pressure,  density,  viscosity,  Reynolds  number,  and  speed  of  sound.4 
This  routine  was  adopted  from  an  IBM  1620  program  for  the  calculation 
of  atmospheric  properties  up  to  2320  kilofeet.  A  polynomial  function 
of  temperature  as  a  function  of  geopotential  altitude  is  used  up  to 
approximately  295  kilofeet.  Below  this  altitude  it  is  assumed  that  air 
has  a  constant  molecular  weight  of  28.  9644.  Geometric  altitude  is 
employed  r^bove  295  kilofeet.  For  these  higher  .  Ititudes  the  moleculai 
weight  is  given  by  the  expression 


Ms  32.  154  +  AKM  [-0.  034513  +  A  KM  '2.  0326  x  10 
+  2.  1032  x  10"8  AKM)] 


(74) 


where  AKM  is  the  altitude  in  kilometers.'  Molecular-scale  temperature 
(°R),  the  rate  of  change  of  temperature  with  geopotential  altitude,  and 
the  initial  pressure  values  (lb-ft  )  are  listed  for  23  altitude  values  in 
the  atmosphere  table. 


The  kinetic  temperature  is  obtained  from  the  molecular-scale 
temperature  through  the  equation 


T  -  Tm  M 
28.9644 


(75) 


where  the  molecular -scale  temperature  can  be  obtained  from  the 
relation 

Tm  (X)  =  D  +  BX  +^Ei  |  X  "  Xi  |  (76) 

i 

The  X^'s  are  altitudes  at  the  end  points  of  straight  line  segments  over 
which  the  variation  of  altitude  and  molecular -scale  temperature  is 
assumed  linear.  The  constants  D,  B,  and  E2  are  functions  of  the  X^ 
and  the  end-point  temperatures  which  are  denoted  as  Yj. 
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(j  =  2.  .  .  23) 

(77) 

D  =  C  +  H2X,  +  H23  X23 

(78) 

c=^ 

(79) 

B  =  -  H,  -  H„ 

(30) 

Ei  =  «j  -  Hj  tl 

(81) 

Once  the  free-stream  conditions  are  known,  the  laminar -turbulent 
transition  coordinate  is  calculated  as  a  function  of  the  free-stream 
conditions  and  using  Equation  (64).  The  axial  viscosity  is  obtained 
from  the  Sutherland  Equation  (62)  fc_  the  laminar  region  and  from 
Equation  (63)  in  the  turbulent  zone.  Constants  A  through  P'  which  appear 
in  the  chemical  production  terms  are  evaluated  on  the  axis  using  p0» 

T0  ,  and  q0  ,  .  .  .  and  at  the  edge  of  the  w~ke  using  pe>  Te,  and  ae»  •  •  •  • 

Density  species  derivatives  in  Equations  (45)  through  (52)  are  then 
calculated  at  the  axis  and  on  the  outer  edge.  From  the  profiles  for 
species  concentrations,  velocity,  and  enthalpy,  a(x,n),  u(x,n),  and 
H(x,  n)  can  be  evaluated.  T(x,n)  is  then  obtained  by  iteration  of  Equation 
(53)  using  the  individual  species  enthalpies  given  in  Equation  (54).  Once 
the  pressure  is  evaluated,  from  Equation  (68),  p  (x,  n)  is  given  by 
Equation  (55). 

At  this  point,  9C  is  calculated  from  Equation  (73).  This  is  only 
done  for  the  initial  coordinate.  Once  6m  is  obtained  from  Equation  (40), 
the  derivatives  of  a0^,  ae  (i),  p,  u0  ,  and  ue  in  Equations  (36),  (38),  (68), 
(34),  and  (37),  respectively,  allow  these  parameters  to  be  calculated  at 
the  next  downstream  coordinate.  Equations  (58)  and  (56)  give  d£n  Te/dx 
and  d£  npe/dx,  and  6  is  obtained  by  integration  of  Equation  (72).  These 
lead  to  d0/dx  Equation  (39).  Ac  and  Hnng  from  Equations  (44)  and 
(43)  allow  dH0  /dx  to  be  determined  using  Equation  (35). 

Finally,  T0  can  be  evaluated  using  the  new  values  of  H0  ,  u0  ,  and 
a0^  ,  and  the  calculations  are  then  repeated  for  an  incremented  value 
of  x.  In  this  calculation,  x  was  incremented  uniformly  as  A  x  =  0.  1. 
Better  results  may  be  obtained  by  using  smaller  increments  for  small 
values  of  x;  however,  the  error  was  not  large  enough  to  justify  a  proce¬ 
dure  of  this  type. 
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Appendix  A 


FORTRAN  PROGRAM 


The  hypersonic  WAKE  program  described  in  the  text  is  given  along 
with  the  required  subroutines.  Comments  appearing  in  the  right-hand 
margin  refer  to  equations  discussed  in  the  text.  Some  important  symbols 
are  defined  in  the  list  of  symbols;  however,  most  of  the  quantities  can  be 
identified  by  reference  to  the  sequence  of  calculation  (Paragraph  6)  and 
the  original  equations. 
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Appendix  B 


SAMPLE  CALCULATIONS 


Sample  cases  involving  a  blunt  body  (4>^=  90  degrees)  and  a  cone 
20  degrees)  were  run  at  150  kilofeet  and  Mco  =  17.  2.  The  calcula¬ 
tions  began  at  x  =  3  and  carried  through  x  =  50.  As  an  example,  the 
figure  compares  the  axial  temperature  as  a  function  ot  x  for  the  blunt 
body  and  the  cone  The  transition  from  laminar  to  turbulent  flow  cjcurs 
sooner  for  the  blunt  body  thereby  inic. ating  a  rapid  decay  in  axial 
temperature. 

Input  data  for  the  cone  are  listed  below.  Velocities  are  in  feet 
per  second.  Temperature  is  in  degrees  Kelvin-  Enthalpies  are  in  feet2 
per  second2.  Altitude,  diameter,  and  x  are  in  feet.  EXOUT  deter¬ 
mines  the  amount  of  output  data.  For  minimum  output,  EXOUT  =.  0.  0 
or  is  omitted.  If  EXOUT  f  0.  0,  intermediate  results  at  each  Ax  =  0.1 
are  printed  as  listed  preceding  scatement  90  in  Appendix  A. 
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